Standard Model Neutrinos as Warm Dark Matter 
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Standard Model neutrinos are not usually considered plausible dark matter candidates because the 
usual treatment of their decoupling in the early universe implies that their mass must be sufficiently 
small to make them "hot" dark matter. In this paper we show that decoupling of Standard Model 
neutrinos in low reheat models may result in neutrino densities very much less than usually assumed, 
and thus their mass may be in the keV range. Standard Model neutrinos may therefore be warm 
dark matter candidates. 
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I. INTRODUCTION 

Cosmological dark matter in the form of neutrinos with 
masses in the eV range is the quintessential example of 
hot dark matter (HDM) (lj . Indeed, for a standard relic 
that decouples from the surrounding thermal bath when 
still relativistic, the current abundance can be easily es- 
timated to be H 
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where gx is the number of degrees of freedom of the 
particle and <7*(Td) is the total effective number of rel- 
ativistic degrees of freedom at the decoupling tempera- 
ture Tjj. Relativistic standard model neutrinos decou- 
ple from chemical equilibrium at a temperature of a few 
MeV, when g^{T£>) ~ 10.75. Since current observations 
indicate that the dark matter density amounts to approx- 
imately fix — 0.3, one usually concludes that the mass 
of standard model neutrinos cannot be larger than about 
30/i 2 eV@§. 

Since neutrinos must be light in order to avoid overdos- 
ing the universe, they were moving at nearly the speed 
of light at redshift z ~ 10 6 when the cosmic horizon first 
encompassed 10 12 M Q , the amount of dark matter con- 
tained in the halo of a large galaxy like the Milky Way. 
This implies that the free streaming of light neutrinos 
destroyed any fluctuations smaller than that of a super- 
cluster (about 10 15 M Q ). Cosmological structure forms in 
a light-neutrino dominated universe in a top-down sce- 
nario, in which superclusters of galaxies form first, with 
galaxies and clusters forming through a process of frag- 
mentation. However, in this scenario galaxies form too 
late and their distribution is much more inhomogeneous 
than observations indicate. An even more serious prob- 
lem is that when normalized to the low amplitude of the 



cosmic microwave background fluctuations detected by 
the COBE satellite |J, the HDM spectrum is only be- 
ginning to reach nonlinearity at the present epoch and 
the free-streaming wavelength cutoff in the amplitude of 
the spectrum should be considerably smaller to form any 
structure by the present. Because of these difficulties, the 
idea of standard model neutrinos as dark matter has been 
gradually abandoned in favor of scenarios where most of 
the mass in the universe is in the form of cold dark matter 
(CDM). 

The purpose of this paper is to show that in low- reheat 
cosmologies, standard model neutrinos may play the role 
of warm dark matter (WDM), avoiding the pitfalls of 
HDM, and may be responsible for the cosmological struc- 
tures we observe. 

The key point of our observation is that the limit to the 
contribution of neutrinos to the present density need not 
constrain neutrino masses to be in the eV region. On the 
contrary, neutrino number densities in low-reheat models 
may be more than an order of magnitude lower than usu- 
ally assumed, and their mass, m v , more than an order 
of magnitude higher, perhaps in the keV range. Fluc- 
tuations corresponding to sufficiently large galaxy halos 
with masses around lO n M0 may survive free streaming 
and standard model neutrinos may act as WDM. Be- 
fore launching into more details, let us see why standard 
model neutrinos may be heavier than usually thought. 



II. ON THE NEUTRINO DENSITY 

The traditional computation of the abundance of stan- 
dard model neutrinos in the early universe is based on 
the simple, but untested, assumption that light neutri- 
nos were in chemical equilibrium at temperatures larger 
than Tjj. In other words, in the standard cosmological 
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30/i 2 eV, it is tacitly assumed that the uni- 



verse had gone through a radiation-dominated phase with 
temperatures larger than about an MeV with active neu- 
trinos in equilibrium. The assumption of an initial con- 
dition of neutrinos in thermal and chemical equilibrium 
in a radiation-dominated universe is then equivalent to 
the hypothesis that the maximum temperature obtained 
during the (last) radiation-dominated era, which will re- 
fer to as the reheating temperature Trh, is much larger 
than the decoupling temperature. The fact that we have 
no physical evidence of the radiation-dominated era well 
before the epoch of nucleosynthesis is a simple, but cru- 
cial, point. It was shown in Ref. || that the reheating 
temperature of order 1 MeV is still compatible with light- 
element nucleosynthesis. A priori, one should consider 
Trh as an unknown quantity that can take any value as 
low as 1 MeV. 

It is usually assumed that the radiation-dominated era 
commences after a period of inflation, and that the cold 
universe at the end of inflation becomes the hot universe 
of the radiation-dominated era in a process known as 
reheating. The reheating process need not be instanta- 
neous. On the contrary, before the radiation-dominated 
phase there may have been a prolonged phase during 
which the energy density of the universe was dominated 
by some component other than radiation. This compo- 
nent is often represented by a coherent oscillating field 
such as the inflaton held, but one could just as easily 
imagine that the universe is dominated by some unsta- 
ble massive particle species. During reheating there is 
a slow formation of a thermal bath of rclativistic par- 
ticles. The temperature of this thermal bath has a pe- 
culiar behavior [2j]. It reaches a maximum temperature 

T max ~ T RH {HjMlJT^jjf'* [Hi is the value of the 
Hubble rate at the beginning of the reheating process) 
and then has a less steep dependence on the scale fac- 
tor a than in the radiation-dominated era, T ~ a~ 3 / 8 . 
During this phase entropy is continuously created (the 
universe is reheating!) and the Hubble rate scales like 
H - [g*{T)/gl /2 (T RH )]{TVTl H M Pl ). 

At a given temperature, the expansion is faster for 
smaller reheat temperatures. When the temperature de- 
creases to Tjih, the universe enters the radiation phase, 
and one recovers the more familiar Hubble law, H ~ 
T 2 /M Pl . 

Let us now assume that the largest temperature of 
the universe during the radiation-dominated phase is 
very small, of the order of a few MeV. Since neutrinos 
have only weak interactions, it is very difficult for the 
thermal scatterings during the reheating stage to gen- 
erate standard model neutrinos through processes like 
e + e~ — > vv and to bring neutrinos into chemical equi- 
librium. Furthermore, decreasing the reheat tempera- 
ture increases the rate of the expansion of the universe, 
making it more and more difficult for the weak inter- 
actions to bring the neutrinos to chemical equilibrium. 
Therefore, if the reheating temperature is small enough, 



standard model neutrinos produced during the reheating 
stage never go into chemical equilibrium. In this case, 
neutrinos are present in the thermal bath at the begin- 
ning of the radiation-dominated phase, but they have a 
number density n v that is smaller than the equilibrium 
number density p). This simple argument shows that 
the present abundance of neutrinos may be much smaller 
than predicted assuming that the largest temperature of 
the radiation-dominated universe was much larger than a 
few MeV. Thus, standard model neutrinos heavier than 
about 30/i 2 eV are perfectly compatible with cosmology. 

The above expectation was confirmed in Ref. jf| , where 
the effective number density of neutrinos was computed 
by solving the corresponding Boltzmann equation ob- 
tained under the assumptions of Maxwell-Boltzmann 
statistics and local thermodynamic equilibrium in the 
calculation of the thermal averaged cross section. This 
amounts to assuming that electrons have an equilibrium 
Boltzmann distribution function and to neglect all the 
Pauli blocking factors. In this paper we have performed 
a more refined and correct computation by solving di- 
rectly the kinetic equations for the neutrino phase-space 
distribution f v (p, t) . 

We will use the numerical code developed in Ref. [0, 
where the collision integrals are analytically reduced to 
two dimensions [0. The kinetic equations for neutrinos 
have the form 
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where the collision integral I co \\ is dominated by two- 
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Here F = / 3 / 4 (l - /i)(l - h) - /i/a(l - /s)(l - h), \M \ 2 
is the square of the weak-interaction amplitude summed 
over spins of all particles except the first one, and S is the 
symmetrization factor which includes (1/2!) for each pair 
of identical particles in the initial and final states and a 
factor of 2 if there are two identical particles in the initial 
state. Finally, there is a summation over all possible sets 
of leptons 2, 3, and 4. Notice that a similar approach 
was taken in Ref. [0 , even though there the interactions 
among neutrinos as well as the electron mass were ne- 
glected and electrons were assumed to have a Boltzmann 
distribution. (These approximations allow the reduction 
of the collision integrals in the kinetic equations to one- 
dimensional integrals, thus simplifying the numerical cal- 
culations.) Furthermore, in Ref. |5|| the emphasis was on 
the impact of very low reheating temperatures on stan- 
dard big-bang nucleosynthesis, and therefore neutrinos 
were taken to be massless. 

We have assumed that reheating is due to the decay 
into light states of a particle 0, which might be the infla- 
ton field, a modulus field, or any unstable particle which 



dominated the energy density of the universe before the 
radiation-dominated phase. The time evolution of the p^ 
energy density is given by 



(4) 



— = -Tptj, - 3Hp^, 



where T is the decay rate of the <j) field, which may 
be expressed in terms of the reheat temperature Trh 
as T = 3H = 3(T| H /Mp i )[87r 3 5 ,(T i? , ff )/90] 1/2 where 
g*(TRn) is the number of relativistic degrees of freedom 
at Trh- Notice that in Refs. the value of 5* in the 

(T-Trh) relation was fixed to 10.75. This is not neces- 
sarily the actual value of 3*(Trjj); we will return to this 
point below. Also, in Ref. H, the definition of Trh is 
r = H, rather than T = 3H. 

We imposed the covariant energy conservation pit) = 
—3H(p + P), where p is the total energy density is 



P = P<f>(t) + 
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and a similar expression holds for the pressure P. The 
initial density of the scalar field p^, is not relevant since 
it just defines the initial time of the evolution. 
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FIG. 1. The relative neutrino densities as functions of time 
for T RH = 1 MeV. 

At the beginning of the evolution, the electromagnetic 
interactions are much faster than the neutrino interac- 
tion rates and therefore, before neutrinos are produced, 
the plasma reaches temperatures higher than the final re- 
heating temperature Trh [iUH- Later, neutrinos Vi of a 
given flavor start being produced by electron-positron an- 
nihilations (e + e~ — > ViVi) and by neutrino- (anti)neutrino 
annihilations {vjVj — > v\V\ with i ^ j). While tau- and 
muon-neutrinos are produced only by neutral current in- 
teractions, electron-neutrino production has a contribu- 
tion from charged current interactions as well. However, 
since the present bound on m Ue is in the eV range, we will 
not be interested in i/ e 's as dark matter. At late times, 



the neutrino distribution reaches some dynamical shape 
which differs from the equilibrium one. We have checked 
that the neutrino abundances follow the same evolution 
curve regardless of their initial abundances and the final 
neutrino distribution is insensitive to the initial condi- 
tions. The time dependence of neutrino density (at late 
times) is shown in Fig. ^. 
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FIG. 2. The relative number density of different neutrino 
species to the equilibrium neutrino number density as a func- 
tion of the reheating temperature. 




FIG. 3. The final distribution functions for v e , and v T 
for Trh = 1 MeV as function of the dimensionless momentum 
Pv/T-y. The dotted line corresponds to the standard model 
neutrino distribution function. 

The final tau- and muon-neutrino number density 
(they are the same) normalized to the neutrino num- 
ber density in equilibrium n v /rJ^ (where — 
2£(3)T !/ 3 /tt 2 ) i s plotted in Fig. § as a function of T RH . 
This ratio becomes smaller than unity for Trh smaller 
than about 8 MeV, signalling a departure from equi- 
librium. The final momentum distribution function for 
Trh = 1 MeV is plotted in Fig. |. For comparison, we 
have plotted the equilibrium neutrino distribution func- 
tion taking into account the fact that the neutrino tem- 
perature is a factor 1.4 times smaller then the photon 
temperature T 7 due to e ± annihilation after neutrino de- 
coupling. 



What is relevant is that the abundance of tau- and 
muon-neutrinos is about a factor of 2.7 x 1CP 2 smaller 
than the standard abundance for reheating temperatures 
around 1 MeV. This in turn implies that the standard 
upper bound of 30/i 2 eV on m„ no longer applies and 
the abundance of tau- and muon-neutrinos for 1 MeV < 
Trh ~ 3 MeV can be expressed as 

^' = V'=£)(^)' W 

where we used a power-law approximation n„/n™ « 
Q.024(Tflij/l MeV) 3 . For larger reheat temperatures, 
3 MeV < T RH < 8 MeV, one needs a slightly 
more complicated expression to fit n v : n v /n^ ~ 
0.44 tan" 1 [(T RH (MeV) - 2.6)/1.17] + 0.43. 
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FIG. 4. The actual value of g, at Trh as a function of Trh- 

As mentioned above, the value of in the (r— Trh) 
relation is actually a function of Trh- The value of 
9*(Trh) provides an indication of the relative abundance 
of neutrinos in the thermal bath. The value of g*(Tmi) 
as a function of Trh is given in Fig. [| and allows to 
relate the decay rate to the reheat temperature, r _1 = 
{MeV/TRH) 2 g*(TR H )- 1/2 1.6 sec ~ 0.56(MeV/7W) 216 
sec. 

Our result in Eq. (^) is slightly smaller than the relic 
neutrino abundance computed in Ref. |w. For instance, 
Fig. 2 of Ref. || shows a value of £l y hr a factor of 1.6 
larger than what is given by Eq. @, for T RH = 2 MeV. 
The calculation in Ref. || makes the approximation of 

(i) neglecting the electron mass in the collision integrals, 

(ii) assuming a Boltzmann distribution for the electrons, 
(Hi) neglecting v-v interactions. We find that the first 
approximation is well justified, but the other two give er- 
rors of up to 30% (depending on the value of Trh), which 
however are opposite in sign and roughly cancel each 
other. The numerical discrepancy is largely explained by 
the fact that the authors of Ref. || fix <?*(!##) — 10.75, 
and therefore their result should be multiplied by a cor- 
rection factor [10.75/<7*(Trjj)] 3 / 4 . For the same reason, 
the limit on Trh derived in Ref. S from nucleosynthesis 



should also be rescaled by a factor [10.75 /g^TRn)] 1 ' 4 - 
With this procedure we find that nucleosynthesis gives 
slightly more stringent lower bounds on Trh ■ Trh > 1.2 
MeV (at 68% CL) and T RH > 0.8 MeV (at 95% CL). 

The result of Eq. (|^) is similar to that found in Ref. || 
(allowing for the differing definition of Trh mentioned 
above). This small difference traces to the assumptions 
made in Ref. || of Maxwell-Boltzmann statistics and 
local thermodynamic equilibrium in the calculation of the 
thermal averaged cross section. 

III. STANDARD MODEL NEUTRINOS AS 
WARM DARK MATTER 

CDM reproduces the observable universe at large 
scales, but appears to be in conflict with observations 
on sub-galactic scales. CDM produces too many dwarf 
galaxies and overdense galactic cores compared to ob- 
servations. The resolution of this difficulty has been 
searched for along different routes. Many suggestions can 
be united under recipe of reducing the power on small 
scales. In the WDM scenario, this reduction occur natu- 
rally via the mechanism of free streaming. The comoving 
smoothing scale can be estimated as the comoving hori- 
zon at matter-radiation equality times the rms velocity 
of dark matter particles at that time 

R sa 0.2 (fi„/i 2 ) 1/3 (keV/m x ) 4/3 Mpc. (7) 

Interestingly, Qx also takes correct value within the 
same range of mx if the reheating temperature is of order 
1 MeV [see Eq. (g)]. This numerical coincidence deserves 
further study and detailed numerical investigation. 
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FIG. 5. The effective chemical potential as a function of 
reheat temperature. 

Encouraging preliminary results were obtained in Ref. 
fLOf . To aid such studies, we parameterize our final dis- 
tribution function of neutrinos as a thermal distribution 
with some effective "temperature" and a "chemical po- 
tential." The difference between the effective tempera- 
ture and the standard-model neutrino temperature is in- 



significant (only a few percent) and the modification to 
the distribution function can be attributed to the chem- 
ical potential, which is shown in Fig. ^|. 

At large Trh the chemical potential approaches the 
standard-model value [i — O.OIT^. At small Trh it can 
be fitted as fj,/T v = 8.8 exp(-0.84T RH /MeV). 

To conclude, the combination of a standard model neu- 
trino of mass of a few keV and a reheat temperature 
of about an MeV will result in neutrinos as candidate 
WDM. Assuming Vl v = 0.3 in one species of neutrino 
and h = 0.65, then R = 0A(T RH /MeV) 4 Mpc. As dis- 
cussed above, any reheat temperature above about an 
MeV is consistent with BBN. 

We want to stress that the hypothesis of neutrinos 
as WDM is not inconsistent at present with the recent 
data on atmospheric and solar neutrino anomalies. At- 
mospheric neutrino data could be accomodated by oscil- 
lations between quasi-degenerate and v T states and in 
such a case solar neutrino observations would require v e 
to oscillate into a sterile state. Alternatively, the solar 
neutrino data might be explained by conversion of v e into 
an active state and the atmospheric neutrino deficit by 
the conversion of into a sterile state. 
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